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Excitons, composite electron-hole quasiparticles, are known to play an important role in optoelec¬ 
tronic phenomena in many semiconducting materials. Recent experiments and theory indicate that 
the band-gap optics of the newly discovered monolayer transition-metal dichalcogenides (TMDs) 
is dominated by tightly bound valley excitons. The strong interaction of excitons with long-range 
electromagnetic fields in these 2D systems can significantly affect their intrinsic properties. Here, we 
develop a semi-classical framework for intrinsic exciton-polaritons in monolayer TMDs that treats 
their dispersion and radiative decay on the same footing and can incorporate effects of the dielectric 
environment. It is demonstrated how both inter- and intra-valley long-range interactions influence 
the dispersion and decay of the polaritonic eigenstates. We also show that exciton-polaritons can 
be efficiently excited via resonance energy transfer from quantum emitters such as quantum dots, 
which may be useful for various applications. 

PACS numbers: 78.20.-e, 73.22.Pr 


Introduction. Monolayer molybdenum disulfide M 0 S 2 
and other group-VI transition-metal dichalcogenides 
(TMDs) are novel two-dimensional (2D) semiconductor 
systems whose electronic and optical properties attract a 
great deal of attention [I-27| . One of their prominently 
discussed features is the opportunity to manipulate the 
valley degree of freedom, including by optical means due 
to opposite-handed circular polarizations of the inter¬ 
band transitions in the two valleys l4|. A growing exper¬ 
imental and theoretical evidence [^tSI I17h 27| indicates 
that the band-gap optical properties of monolayer TMDs 
are dominated by relatively tightly bound electron-hole 
pairs, excitons, with binding energies substantially larger 
than in the majority of conventional inorganic semicon¬ 
ductor quantum wells [ 2 ^, [ 2 ^. The corresponding 2D 
exciton physics in TMDs may therefore reflect generally 
stronger interactions of excitons with macroscopic elec¬ 
tric fields and light. 

An important fundamental issue is the nature of in¬ 
trinsic excitonic eigenstates and their energy-momentum 
(dispersion) relationships. It was shown recently [ 2 ^ that 
the long-range exchange Coulomb interaction mixes in¬ 
dividual valley excitons to establish excitons with the 
longitudinal (L) and transverse (T) polarizations as nor¬ 
mal system modes, similarly to quantum well excitons 
[301 l3Ij . The resulting exciton spectrum, as a function 
of its center-of-mass in-plane wave-vector k = (fcx, ky) = 
k (cos 9, sin 9), was found to exhibit a specific Dirac-cone- 
like behavior at low momenta hk, in particular within the 
light cone, k < q = ui/c. As we perform in this Letter 
a semi-classical analysis of the interaction of monolayer 
excitons with long-range electromagnetic fields, it will 
be shown however that (1) the dispersion of L- and T- 
excitons is affected via both inter- and intra-valley pro¬ 
cesses leading to the overall generic behavior character¬ 
istic of the 2D excitons 33-33| without a Dirac-cone fea¬ 


ture; (2) Moreover, also in a general fashion, the intrin¬ 
sic behavior within and in the vicinity of the light cone 
is that of the 2D exciton-polaritons as determined by 
the full electromagnetic (rather than just electrostatic) 
mixing of valley excitons taking account of the retarda¬ 
tion effects [ 29 I, [SSj. (In the electrostatic limit the ob¬ 
tained dispersion reproduces the exciton spectrum that 
we derive from inter- and intra-valley exchange Coulomb 
interactions). General features of the intrinsic exciton- 
polaritons are clearly accentuated in the macroscopic 
electrodynamics framework, which allows for a straight¬ 
forward generalization of the analysis of the polariton 
dispersion and radiative decay in free-standing monolayer 
TMDs to monolayers at the interface between different 
media. We also use that framework to illustrate the pos¬ 
sibility of very efficient direct excitation of monolayer po- 
laritons by energy transfer from proximal electric-dipole 
emitters such as quantum dots, which may be useful for 
various applications. 

Exciton-polaritons from effective dipole-dipole interac¬ 
tions. In a common basic description of optically active 
lowest energy excitons in monolayer TMDs, they arise 
as a result of the direct Coulomb attraction between an 
electron and a hole in the same valley. Due to the strong 
spin-orbit coupling in TMDs, the electron-hole spin com¬ 
position is associated with the valley and therefore 
will be omitted. In the valley-based picture there are 
then two exciton species (index a = 1,2) correspond¬ 
ing to two different valleys, both having the same energy 
dispersion 


Eo(k) = hoJo(k) = hojo -b fi^fc^/2Mo, 


( 1 ) 


where Hujo is typically somewhat below 2 eV and the 
exciton mass Mq is close to the free electron’s me 
SEIillil. The parabolic kinetic energy in Eq. ([T]) 


with Mo being the sum of electron and hole effective 
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masses signifies the fact that the exciton propagation in 
space is enabled by the simultaneous motion of the elec¬ 
tron and hole. The inclusion of the standard exchange 
Coulomb interaction results in another mechanism of the 
exciton propagation and concomitant modifications of its 
dispersion, which can be seen as the electric field of the 
exciton annihilated at one spatial point creating the ex¬ 
citon at another point. Such a long-range process is well- 
known to correspond to the electrostatic dipole-dipole 
coupling, which is the major exciton transfer mecha¬ 
nism in molecular systems (33|. This physically attrac¬ 
tive semiclassical real-space picture can then be readily 
extended to the dipole-dipole coupling mediated by the 
full electromagnetic interactions with the retardation ef¬ 
fects in place. The resulting intrinsic excitations, exciton- 
polaritons, would thus take into account the exciton-light 
interaction absent in the picture of Coulomb excitons de¬ 
rived with electrostatic interactions alone [ 2 ^, [H, [13 ■ 
The electromagnetic interactions in monolayer TMDs 
may be classified as intra-valley {a = /3) and inter-valley 
(a ^ P) couplings, in some analogy with bipartite lat¬ 
tices and molecular systems with two molecules per unit 
cell [13 ■ Correspondingly, the exciton self-energy correc¬ 
tion due to the interactions is a 2 x 2 matrix S. In the 
continuum description of 2D Wannier-Mott excitons, it 
can be written in the form of 

j \ (2) 


as a function of energy E = hw = hcq and wave-vector k 
variables. Here Vap{r) = |V’(0)P Ki, 3 (r) is the 2D energy 
density determined by the probability |^/’(0)p = 8/(7ra^) 
of finding the electron and hole of the exciton at the same 
spatial point [j^ l and the interaction matrix elements 
V'Q; 3 (r) dependent on the relative two-dimensional posi¬ 
tion r = rf. From Ab initio calculations [s, IT, 21, 
the exciton Bohr radius in monolayer TMDs is estimated 
as as ^ 1 nm. For the free-standing monolayers in vac¬ 
uum the long-range part of the interaction 


piqr , 

V^f3ir) = — {^ 


%q 

r 




d/o.d°^-3(f.d/°)(f.d°^)l (3) 


as arising from the full electromagnetic dipole-dipole cou¬ 
pling (in Gaussian units) [13. The familiar electro¬ 
static limit corresponds to the light vacuum wave number 
q = w/c —>■ 0 (speed of light c —>■ 00 ). The interacting 
a and /3 species in Eq. ([3]) are represented by the cor¬ 
responding interband dipole transition matrix elements 
for creation, d^°, and annihilation, d^^, of electron-hole 
pairs. In the electrostatic limit, the self-energy correc¬ 
tions © would be real-valued and functions of wave vec¬ 
tor k only. With the retarded electric fields, however, 
Sq ,/3 are functions of both energy and wave-vector vari¬ 



FIG. 1: Functional behavior of the (a) real and (b) imagi¬ 
nary parts of the intrinsic exciton-polariton self-energy as per 
Eqs. ((6| and (IIOII . Red color lines are used for the longitudi¬ 
nal (L) and blue for the transverse (T) excitons. Shown with 
dashed lines are the results for a free-standing monolayer in 
vacuum, and solid lines are for the monolayer on a glass-like 
substrate, see text. 


ables and generally complex-valued; they would there¬ 
fore determine both the exciton dispersion and the decay 
width (decay into photons) in a self-consistent calcula¬ 
tion. Calculations with the retarded fields are equivalent 
[33 to calculations of exciton-polaritons as arising from 
the interaction of Coulomb excitons with transverse pho- 

2i[33. 


tons 


The in-plane interband dipole transition matrix ele¬ 
ments idy°‘^'^ have a common form for both 

valleys but with the opposite handedness of their circular 
polarization: 


4 ^^ = df'Tp = 4 ^^ = “ 4 ^^ = 


( 4 ) 


where we chose 2? as a real positive quantity. From a two- 
band monolayer TMD model [3, for instance, the dipole 
transition moment would be 2? = eat/ Eg, where e is the 
fundamental charge, a the lattice structure constant (~ 
3.19 A), t the nearest neighbor hopping integral (~ 1.10 
eV) and Eg the energy gap between the conduction and 
valence bands that the excitons come from. A calculation 


based on the three-band model |23| gives the same form, 
but with a slightly different value of t. 

The difference between the transition dipole moments 
o in the two valleys translates into the difference of 
inter- and intra-valley couplings ([3|). A direct calcula¬ 
tion of integrals in Eq. ([2]) then leads to the following 
k-dependence of the self-energy matrix [33 : 


S(A,k) = 


Jo exp(—2i0) Ji 

exp(2zd) Ji Jq 


( 5 ) 


where the intra-valley component Jq = —irii^/q^ — 4 - 

q^/\/q^ — hP) and the inter-valley component Ji = 
i'qk'^ / \/q^ — feature the same magnitude scale rj = 
27r|V'(0)p2?^. The functional form of these components 
illustrates the great qualitative distinction resulting from 
the retarded interactions: the obtained corrections Jq 
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and Ji are purely imaginary above the light line, k < q, 
but become purely real {\/q^ — k? —> i\/k^ — q^) below 
the light line, k > q. As is known [s^, [s^, this signifies 
the impossibility for the intrinsic 2D exciton-polariton 
with k > q to decay into a photon due to the momen¬ 
tum conservation. Exciton-polaritons with k < q, on the 
other hand, exhibit the radiative width due to such a 
decay. Moreover this width is grea tly enhanced in com¬ 
parison with localized emitters |33l[ . It is also instructive 
to look at the electrostatic limit {q —>■ 0 ) of these correc¬ 
tions: the inter-valley component then becomes Ji = iqk, 
precisely the result derived in Ref. [ 2 ^ . Importantly, the 
intra-valley component Jq becomes equal to Ji in this 
limit, which is also confirmed in our independent calcu¬ 
lations with the exchange Coulomb interaction. 

Because of the inter-valley coupling Ji, the valley ex¬ 
citations are clearly not the eigenstates of the system. 
Instead, the eigenstates are their linear combinations 

'k± = 2 “^/^ ^ ±exp( 2 id) ) diagonalize matrix S 

and yield the corresponding self-energy corrections as 


= Jo ± Ji = —2ir] X 


{L), 

(T). 


( 6 ) 



FIG. 2: Acceleration of the decay rate F of a randomly ori¬ 
ented electric dipole emitter such as qnantnm dots in the 
vicinity of the monolayer on a glass-like substrate as com¬ 
pared to the spontaneous decay rate Fq in vacuum, Eq. (HU. 
The data are shown as a function of the emitter frequency 
Lj and for different distances h to the interface: 5 nm (black 
lines), 10 nm (red) and 20 nm (blue). Solid lines display re¬ 
sults for the nearly dissipationless excitons (7 = 0.1 meV), 
dashed lines are for the dissipation parameter 7 = 25 meV. 
These illustrative calculations were done with the model pa¬ 
rameters dvryo = 2 nm, = 0.25 eV^-nm in Eq. m and 

Ajq = 1-9 eV, Mo = me in Eq. ([T|- 


These eigenstates have, respectively, longitudinal, for 
and transverse, for T* , polarizations with respect to 
polariton wave-vector k 22 |: for the propagation along 
the x-axis, e.g., Eq. (|1]) shows that 4'+ combination corre¬ 
sponds to the transition moment oc {V, +iT>) -\- {T>, —iD) 
along X while 'F_ to the moment oc (2?, -\-iD) — {V, —iV) 
along y. The functional forms of self-energy corrections 
([ 6 |) for the free-standing monolayer are shown in Eig.[T]by 
dashed lines. It is clear that, in the electrostatic limit, 
only the L-excitons would acquire the additional term 
S-|_ = 2rik in their dispersion due to long-range exchange 
interactions, while the dispersion of T-excitons would re¬ 
main unchanged (E_ = 0) as Eq. ([T|)- There is no Dirac- 
cone-like behavior even in the electrostatic limit. In the 
valley-centric basis these conclusions thus follow from the 
simultaneous account of both inter- and intra-valley pro¬ 
cesses as opposed to the inter-valley coupling alone [ 22 1 . 

Exciton-polaritons from macroscopic Maxwell equa¬ 
tions. The obtained results for polaritonic eigenstates are 
in agreement with the picture known for quantum wells 
and reflect the fact that the opposite-handedness in¬ 
plane susceptibilities Xi and X 2 of individual valleys 
just add up in the overall isotropic electrodynamic re¬ 
sponse of the monolayer. The latter is then characterized 
by the scalar susceptibility y defining the monolayer 2D 
current density j = — dTrfwyE induced by the in-plane 
electric field E. It is this current density that enters the 
boundary conditions for the macroscopic Maxwell equa¬ 
tions determining the effects of long-range fields on the 
system excitations 2^, 331. Eor well-separated excitonic 
states (HI), e.g., the 2D scalar susceptibility acquires a 


familiar single-oscillator form 

X(w,fc) = xo + ^/(wo(^) - - 2170 ;) , (7) 


where xo is the background term due to higher-frequency 
transitions and 7 the phenomenological dissipation pa¬ 
rameter. A many-oscillator form could be used to include 
even more specifics for different TMD monolayers. 

The versatile macroscopic framework can be easily ap¬ 
plied to various environments. Here we exemplify this by 
considering the model of an infinitesimally thin planar 
monolayer between two non-magnetic media with dielec¬ 
tric constants ei and £ 2 - One can solve for the eigen- 
frequencies of the system directly. Alternatively, as also 
useful for other problems, one looks at the poles of the 
reflection coefficients for electromagnetic waves in our 
sandwich configuration. The 3D setup involves not only 
the in-plane components of the fields and wave vectors 
but also their z-components perpendicular to the planar 
interface. With the boundary conditions of the polariz¬ 
able interface monolayer , one easily derives the reflec¬ 
tion coefficient amplitudes for p- and s-polarized waves 
as (w, A:)-dependent 



and 


£L 

kzi 



-1 


( 8 ) 


= {kzi - kz2 + ^Mq^x) [kzi + k^2 - ^Mq^x) ^ ■ 

(9) 
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The normal components k^i = isiq^ — of the 

waves in the respective media appear, as usual [s^, re¬ 
lated to the in-plane wave number k for given frequency 
ui = cq. The in-plane component of the electric field in 
a p-polarized wave is along the in-plane wave vector k 
whereas in an s-polarized wave they are perpendicular 
to each other. Hence, the poles of Eq. ([ 8 ]) determine the 
dispersion and decay of the L-polaritons while poles of 
Eq. dH) those of the T-polaritons. Writing the pole equa¬ 
tions in the form of 


-= —iiri X 

X 


{ei/k^i+£2/kz2) \ {L), 
q'^/{kzi + k^2), (T), 


( 10 ) 


one recognizes that the functional dependence in the 
r. h. s. of Eq. (ITUl) reduces to that in Eq. dH]) for the free¬ 
standing monolayer. It is also clear that with the neg¬ 
ligible dissipation and screening due to higher-frequency 
transitions, —1/x — 2ujo{uj — uiQ{k))/A becomes just pro¬ 
portional to the self-energy corrections in the vicinity 
of the resonance, uj loq. The derivation with the ef¬ 
fective dipole-dipole interactions thus fully agrees with 
the macroscopic electrodynamics result, and we obtain 
A = 2u}[)'q/'Kh by comparing Eqs. dH) and (HOI) . Note that 
the screening by the surrounding can substantially af¬ 
fect the exciton radius and binding energy reducing thus 
“the strength” A of the resonance O; the calculations of 
exciton binding are outside of our scope. 

Figure [I] illustrates the differences in the behavior of 
self-energy corrections for two systems: the symmetric 
sandwich with £i = £2 = 1 (a free standing mono- 
layer) and and the asymmetric sandwich with £i = 1 
and £2 = 2.25 (glass-like substrate), as follows from 
Eq.dini). It is transparent that in the symmetric case, the 
T-polariton branch exhibits splitting at the light cone, 
fc —>■ q, as consistent with the divergence of the radia¬ 
tive decay rate dS]) at fc —>• q — 0. Qualitatively differ¬ 
ently for the asymmetric sandwich, there is no divergence 
in the decay rate and the dispersion of the T-polariton 
branch is continuous. Figure HKb) also shows the exten¬ 
sion of the region of the radiative decay: polaritons with 
q < k < y^q can now decay into photons that propagate 
only in the substrate. Beyond that region, however, the 
intrinsic (without scattering effects) polaritons become 
non-radiative and, conversely, cannot be directly excited 
by plane-wave photons incident on the monolayer from 
the infinity, the situation similar to the excitation of sur¬ 
face waves like surface plasmons 32, ^.1^. 

Application to energy transfer. The direct excitation 
of non-radiative modes is however known to be possi¬ 
ble by the near electromagnetic fields. This may be ac¬ 
complished, e.g., with special optical geometry setups 
[H, in, m or via resonance energy transfer (FT) from 
proximal quantum emitters such as molecules or quan¬ 
tum dots 38 -4^. We note that ET from OD emitters to 
(quasi) 2D absorbers has been of increasing theoretical 
and experimental interest, including organic and inor¬ 


ganic semiconductors graphene and M 0 S 2 [431 
(and references therein). ET provides a new decay chan¬ 
nel and is thus manifested in the acceleration of the ob¬ 
served emitter’s photoluminescence decay. 

Using the macroscopic electrodynamics formalism de¬ 
veloped for such applications 38, li^, the decay rate T 
of the randomly-oriented electric-dipole emitter in the 
medium with dielectric constant £1 at distance h from 
the planar interface is derived [3~' 


as 


r/ri = l -h {l/2kl)Re / {kdk/k^i)e 
Jo 

X (^{2k^ - kj)r^p'> + kfr^^'^^ , 


2ikzih 


( 11 ) 


where Ti is the spontaneous decay rate in the uniform 
medium and k\ = £iq^. The integration in Eq. m 
over all values of the in-plane wave-vectors k signifies 
that the near-fields of the emitter are included. Figure 
H] illustrates the results following from Eq. ([TT]) with the 
reflection coefficients ([ 8 ]) and (|9]) for the considered ge¬ 
ometry. The analysis shows that a very large effect ob¬ 
served in Figure [2] is predominantly due to the excitation 
of L-polaritons (poles of similar to a very efficient 
excitation of surface plasmons by ET [H, ; the disper¬ 

sion of the peak with distance h is also clearly seen. For 
the computational illustration here we used a more mod¬ 
erate value of the resonance strength A in Eq. and 
its larger values would lead to even faster ET into the 
monolayer. Figure [2] also shows how the narrower polari- 
ton peaks are spread by exciton damping ( 7 ) processes. 
Excitation of TMD monolayer polaritons via high effi¬ 
ciency ET from neighboring emitters mi ght be attractive 
for various optoelectronic applications |4^ . Its ex¬ 
perimental studies for different emitter frequencies and 
distances could be used for quantification of polariton 
properties. 

Discussion. While our discussion in this paper has 
been focussed on the intrinsic polaritons, it is recognized 
that various scattering processes due to polariton interac¬ 
tions with phonons and defects can significantly affect the 
observable optical properties of realistic TMDs, and this 
appears as an important topic for future studies. Just as 
with the quantum well polaritons, e.g., their thermaliza- 
tion may lead to a specific temperature dependence of the 
radiative lifetime. With the substantial contribution to 
polariton dispersion in TMDs from the long-range inter¬ 
actions we discussed, that temperature dependence may 
deviate from the dependence resulting from the purely 
parabolic exciton dispersion ( 4 ^ . 

Finally, we estimate the experimentally relevant pa¬ 
rameters for M 0 S 2 . Following Ref. [l^, rj would be esti¬ 
mated as ^ 0.75 eV-A, resulting in substantial contribu¬ 
tions to the exciton dispersion. With this estimate, the 
energy unit in Fig. |l] 2 ? 7 q ^ 1.5 meV for q = w/c ^ 0.01 
nm“^. Likewise, one would obtain iirh^A ~ 1.1 eV^-nm 
for the unscreened value of the parameter A in Eq. ©• 
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